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, Considering quarks as the quasi-particles of the model Hamiltonian with four-fermion interaction 

I we study response to the process of filling up the Fermi sphere with quarks. 

(N 

^ ' In the present paper we analyse the process of filhng up the Fermi sphere with quarks modelhng 
the corresponding Sletter determinant in the direct way. The quarks are considered as the quasipar- 
ticles of the model Hamiltonian with four-fermion interaction. In Ref. [1] the ground state of this 
model Hamiltonian was studied in detail and the singularity of mean energy functional depending 
r-| . on current quark mass was found out. In the course of this study we need to trace back the al- 
terations of respective dressing transformations. It will be seen, from what follows, this problem is 
Ph' fairly difficult to treat with well developed machinery of the Green functions [2J in which a system 
^ ■ is described by the Hamiltonian with additional term controlling the particle number fi q'-f^q {fx is 
the chemical potential of bare quarks). However, the quark chemical potential conception runs with 
^ . some uncertainties in interpreting a baryon chemical potential [3] and quantitatively produces the 
^ ! arbitrariness of order 50 - 100 MeV for the phase diagram. Furthermore, each kind of quarks requires 
' its particular chemical potential. Here we also consider the impact of filling up the Fermi sphere on 
the process of quasi-particle formation. Actually, the approach proposed should be appended by an 
analysis of the bound states infiuence but here we formulate this problem only doing hope to get a 
^ . justification of such a position while completing our consideration. 

00 ! We start with briefiy reminding how the quasiparticle concept appears in this approach. In Ref. 
^ I [1] the ground state of Hamiltonian with four-fermion interaction 
> ■ 

>< : n=-q V + im)q-q t^j, q J dy q' t\ q' {A'^AD (x - y) , (1) 

taken in the form of a product of two colored currents located in the spatial points x and y connected 
by a formfactor was investigated and here q = q{x), q = q{x), q' = q{y), q' = q{y) — are the (anti- 
) quark operators. 



qai[x) 



dp 1 



a(f>, s, c) Uai{p, s, c) e'P^ + s, c) v^iip, s, c) e 'P^] , (2) 



where pi = — — m^, i — is the colour index, a is the spinor index in the coordinate space, a"*", a and 
6+, b are the creation and annihilation operators of quarks and anti-quarks, a |0) = 0, 6 |0) = and 
|0) is the vacuum state of free Hamiltonian. Everywhere that a summation over indices s and c is 
meant, the index s describes two spin polarization of the quark and index c plays the similar role for 
the colour. t° = A°/2 are the generators of colour gauge group SU{Nc) and m is the current quark 
mass. The Hamiltonian density is considered in the Euclidean space and 7^ denotes the Hermitian 
Dirac matrices, /x, z/ = 1, 2, 3, 4. {A'^A'l){x—y) stands for the formfactor and its form will be discussed 
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later. The effective Hamiltonian ([T]) arises as a result of averaging the system of quarks influenced 
by intensive stochastic gluon fleld A", see Ref. [1] (for example, by an (anti-)instanton ensemble). 
Based on the form of induced four-fermion interaction the ground state of the system was searched 
as the Bogolyubov trial function composed by the quark-anti-quark pairs with opposite momenta 
and with vacuum quantum numbers i.e. 



W)=T\0), 

T = Up^s exp{ (p [ a+(p, s) b^{-p, s) + a(p, s) b{-p, s) ] } . 



(3) 



In this formula and below, in order to simplify the notations, we refer to only one complicated index 
which means both polarization indices (spin and colour). The parameter (p{p) describing the pairing 
strength is determined by the minimum of mean energy 



E = {(j\H\o) . 



(4) 



Introducing the dressing transformation T we deflne the creation and annihilation operators of quasi- 
particles as A = T a T^^, = T b^T^^. Now the quark field operators are presented in the following 



form 



q{x) 
q{x) 



dp 



dp 1 



A{p, s) U{p, s) e'P'' + 5+(p, s) V{p, s) e^P"" 
A+{p, s) U{p, s) e-'P"" + B{p, s) V{p, s) e^^^ 



and the spinors U and V are given as 

U{p, s) = cos{{p) u{p, s) — sin(y9) v{—p, s 
V{p, s) = sm{ip) u{—p, s) + cos{{p) v{p, s 



(5) 



where U{p, s) = U~^{p, s) 74, V{p, s) = V~^{p, s) 74 are the Dirac conjugated spinors. 

Then the problem of our interest can be formulated in the following way. We need to construct 
the state filled in by quasi-particles (the Sletter determinant) 



\N)= n A''iP;S)\a) 

\P\<Pf;S 



(6) 



which possesses the minimal mean energy over state |A^). The polarization indices run here over all 
permissible values. The momenta and polarizations of quasi-particles filling in the Fermi sphere are 
marked by the capital letters. The free Hamiltonian written in the quasi-particle operators looks like 
in the following 



Hq = — J dx q{x) (z7V + im) q{x) 



(2vr) 



3 \pi\ [cose A^{p; s)A{p; s) + sin 6' A+(-p; s)B^{p; s) + 
+ sin 6* B{-p; s)A{p; s) - cos 9 B{p; s)B^{p; s) 



(7) 



There are two diagonal matrix elements of free Hamiltonian and now we are going to consider one 
of them 

{N\ Bip;s) B+ip;s) \N) ~ {a\AiP;S) B{p-s) B+ip;s) A+(P;5)|a) . 

While presenting the matrix element we demonstrate in the right-hand part only one partial con- 
tribution constructed by operator A of some detached sort. The others as it can easily be seen 



2 



give a sum (integral) over states filling the Fermi sphere in. Taking into account the normalization 
condition which we assess as {(t\A A'^\a) = 1 (for operators with coinciding arguments A, A~^), the 
matrix element leads to the vacuum contribution similar to one known from Ref. [1] 

- / ^ {N\ \p,\ cose B{p-s) B+{p-s) \N) = - j \p,\ cose . 

We should also remember that B\a) = 0, A\a) = 0. The partial contribution to the second matrix 
element equals to 

{a\A{P;S) A+{p;s) A{p; s) A+{P;S)\a) = (27r)3 S{p - P) S,s {^\A{P; S) A+{p;s)\a) . (8) 
Having filled the state this contribution occurs as 

/ (1^ """"^ ^^(p;^) ^(p;^) 1^) = b4| cos^ . 

There are other diagonal matrix elements in the interaction Hamiltonian q t°"^^ q q' t^'ju q' which 
we are going to designate as 1), 2), a), 7), S). For example, 1) is the following matrix element 
{N\B B+ B' B'+\N), 2) corresponds to {N\B A A'+ B'+\N), a) looks like {N\A A+ A' A'+\N), 
13) is presented by {N\A A+ B' B'+\N), 7) corresponds to {N\A+ B+ B' A'\N) and 5) looks like 
{N\B 5'+ A'+ A'\N). The contribution 1) to the interaction term (A^| q ^7^, q q' t^7^ q' \N) leads 
to the following spinor form 

Bearing in mind the completeness property of spinor basis there appears the unit matrix while 
summing up all colour polarization, for example, X^c ^(c)V^j(c) = and, therefore, the contribution 
1) as was mentioned in Ref. [I] turns out to be zero. The partial contribution 2) looks like 

{a\A{P-S) B{p-s) A{q-t) A^{p'-s')B+{q'-t') A+{P- S)\a) = 

= (27r)6 [5{q - p')5ts'{o\A{P- S)A-^iP; S)\a) - 5{q - P)5ts{(T\A{P; 5)A+(p'; s')\a)] 5{p - q')5st>. 

One can obtain from this expression that it contributes to the matrix element {N\qt°''~^^ q q' t^ju q' \N) 
(if the normalization condition of one single state is taken into account) with the following result 



-V.,{p; s) U^j{P; S)U,,{P; S) tl^^s Vsi{p; s) + 



+V^i{p;s) t^^.7^^ Up,iq;t)U^kiq;t) tliYys Vsiip; s) . 

The calculation of matrix elements (3), 7) and 6) demonstrates that they have a similar structure 
but their contributions cancel each other. 

The matrix element a) deserves the special discussion. Unlike the aforementioned examples the 
major partial contribution here should be searched for the pair of quasiparticles 

{a\A{Q;T)AiP;S) A+(p;s) Aiq;t) A-^{p'; s')A{q';t') A+(P; 5)A+(Q; T)|a) . 

When the momenta Q and P coincide, as known, (see, for example, Ref. ^) the next term in the 
1/V decomposition {V is the volume occupied by the system) appears. But if we are interested 
in knowing how one quasiparticle modifies the dressing transformation it is necessary to consider 
its matrix element with |1) = A'^{P; S)\o'). Omitting the intermediate calculations we present the 
contribution of the scattering term a) as 

-Ua^iQ■, T) tl^r^p Up,{P- S)U,,iP; S) tla^s UsiiQ; T) + 

+U^UP; S) Up,{P- S)U,,{Q; T) tla^s UsiiQ; T) + 

+UM T) U^Ap; s)U,u{p- s) tla^s UsiiQ; T) + 

_ _ 9 

+U_M; T) t^^Y^^ Up,{Q- T)U,,{P; S) tla^s Usi{P; S) + 

-U^P; S) t» 7::^ ^ft-(Q; Tp.kiQ; T) t^a^s Usi{P; S) + 

+UUP- S) t^^^^p Up,{p- s)U,k{p; s) tla^s UsiiP; S) , 



The second and fourth terms give zero contribution for the same reasons which were given for the 
contribution 1). The first and fifth terms, up to the re-designations, are equal to each other. The 
third and sixth terms presents the contribution of two states with momentum P and Q. It is not 
difficuh to understand that for the ensemble we get simply the integral over the Fermi sphere (see 
below). 

When one quasi-particle is considered the matrix element unlike (l9l) is proportional to 

U^^iP; S) t%^% Up.ip; s)U,kip; s) tla^s UsiiP; S) . (10) 
The polarization matrices of quasi-particles were considered in Ref. [1] and have the forms 



VV = P474 + cos(^^) (p7 — im) — 
UU = P474 + cos(6') (p7 + im) + 



a + a 
2im 

a + a 
2im 



sin(^^) (p^ — im p7) 
sin(^) (p^ + im pj) 



where the angle 6 = 2(f, (the summation over polarization was done and the matrix is diagonal in 
the colour indices). As a result for the traces which we are interested in we receive 



Tr {UU 7^ U'U''yu t") J 



I 2 r I „ 01 + a p \ ( . a + a q 

X < ^4^4 gau ± "m Onp cos 9 5- sin 6' cos 9 5- sin 9 

' ' 2 m I \ 2 m 



Ill ) 



a + a , \ I ^, a + a' , ^, , 
cos 9 H — sm 9 q^pi + cos 9 H sm 9 p^qi 



cost 



« + « ■ n\ I /)/ , 0^+0^' . n> 

smy cost/ H smf 



+ 



where p = \p\, q = \q\ and 9' = 9{q). The coefficient a is defined up to the phase factor as a.a= rr? jjp' 
which is quite natural because a spinor has such a phase degree of freedom. The analysis performed 
in Ref. [T] demonstrates that the minimal value of mean free energy is reached when this coefficient 
becomes the real number a = j^- -^^^ clarity we choose the positive sign and specify the form of 

correlation function which is motivated by (anti-)instanton ensemble (for details see Ref. [1]), 



{AI.ADix -y)= S'^" G [I{x - y) ^ " J,A^ " v)] 



(12) 



Here the second term is spanned onto the vector of relative distance. 

Now we have to collect together all the results obtained for the one quasi-particle. For the matrix 
element of interaction Hamiltonian we have 



(l|gt"7,.gg lu q\l) 



2 4|p4||g4| 



xTr 



-V{p)V{p)^,U{P)U{P)^, + V{p)V{p)^,U{q)U{q)^, + U{P)U{P)^,U{p)U{p)^, 



(13) 



The polarization indices are omitted here. It is interesting to notice that in the first and third 
terms the factors spanned onto the tensors g^y, (^4/i^ii/ + ^^s^v^i^), {5i^5jy — 5ij5^i, + Sii,6j^), (see Eq. 
fllip ) cancel each other. They survive only in the second term which describes the pure vacuum 
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contribution (considered in Ref. [T], see also discussion below). Let us define the partial energy 
density per one quark degree of freedom, as 

£ 



w 



2iV. 



£ = E/V 



where E is the total energy of ensemble. Collecting all contributions together we have for the state 
with the single quasi-particle 



■^^1 = Ip^I cos9 + 2G J 



iP^ (^^ ^ - - cos ^"j f sin e'-- cos ^^'^j (/ - J/4) + w^^c , (14) 

{2ny \p4\\q4\ V P J \ (1 J 



where I = I{p + q), Jij = Jijijp + q), J = Y^i=i Ju- The term w^ac describes the vacuum contribution 

IG 

and is shown below. It is convenient to pick out the colour index G = The energy of the state 
develops a minimal value if the condition 

dwi 



dd 







is satisfied. Neglecting the modifications in the vacuum contribution w^ac we have 

- |P4| sin0 + 2G ^ l^sin^-^cos^^l^;^ ^ fsin^'--cos0') (/ - J/4) = , (15) 



(27r)3 \q^\ 



which with the precision up to the terms gfj,^, {Si^Sju — + 6iuSjfj), is equivalent to the condition 

of minimum of mean vacuum energy ^^^"'^ = 0, see Ref. [1], i.e. in the situation of a single quasi- 
particle the dressing transformation in this approximation does not change compared to the vacuum 
one. 




200 400 

Pp MeV 

Figure 1: The dynamical quark mass (|Mq|) as a function of the Fermi momentum for the NJL 
model. The solid line corresponds to the current quark mass m = 5.5 MeV. The dashed one shows 
the behaviour in the chiral limit. 



For the matrix element of the interaction Hamiltonian for the ensemble of quasi- particles we 
have 



(iV| gt»7^gg'r 7,g' \N) 



2 4|p4||g4| 
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xTr [-V{p)V{p) 7^ U{P)U{P) 7. + V{p)V{p) 7^. U{q)U{q) 7.+ (16) 
+U{P)U{P) 7^ f/(p)I7(p) 7. + 2 U{P)U{P) 7, U{Q)U{Q) 7.] , 

(see Eq. (1131) to compare). It is pertinent to mention the coefficient 2 in front of the last term. 
We present the matrix element implying some ordering of the momenta, for instance, |Q| < \P\. 
But in the expressions which are averaged over the state |A^) it is convenient to present the result 
by integrating over the whole Fermi sphere, without taking into account the ordering. It is evident 
then that we should take the contribution in factor two smaller in the interaction term, see the 
formula below. In the last term similar to the 'vacuum' matrix element (the second term there) the 
contributions spanned on the tensors g^^, {5i^5jy — 5ij5^u + 5iy5j^) survive. For simplicity we consider 
in this paper only the situation when the second correlator equals zero J^^ = 0, besides we neglect all 
the distinctions provoked by the tensors Qf^,^, {5ifj,5jy — 5ij5^u + 5iy5j^). Collecting all the contributions 
together for the mean partial energy up to an unessential constant we have 



, Pf dv 

{N\w\N)= I ^ |p4|cos^ + 



fPp dp p m \ [Pf dq q 

-G / TTT ■; Smb* cosy / 7 TTT -; — 

{2tiY \pA \ pjJ (27r)3 \q. 



+2Grj%J^ sin^-^cos^W^^ sin^^'-^cos^' 

(17) 




f dp . ..^ „, ^ f dp P . n ^ n\ f dq q . nl ^ nl\ T 

+ / TTT P4 (1 — COSC'I — G / 7 1 r Smt^ COSO' / -; r Sm COS 1 . 

J (27r)3 l^'l^ > J (27r)3 \p^\ \ p J J {27if \q4 \ q J 

To obtain this we perform a natural regularization subtracting the contribution of free Hamiltonian 
Hq. As a result there appears an unit in the bracket containing — cos^ in the last line. Performing 
the transposition of integration and changing the variables as p — q the interaction term can be 
rewritten in the following conventional form 2 J^'^dp J dq — J^^dp J^^ dq — J dp J dq = — Jp^dp Jppdq. 
Finally we obtain the following expression for the partial energy 



, dp P ( . ^ m \ f dq q ( . ^, m \ 
G / ——IT -; — r sm cos / ——^ — r sm cos ] I . 



Ipp (27r)3 \p^\ \ p ) Jpp (27r)3 \q^\ \ q 

It has quite interesting interpretation. Compared to the vacuum mean energ}0, see. [1], 

dp I .^ ^ f dp P (. n ^ n\ f dq q ( . m / 
1^4 1(1 — COS 0) — G - — -o- -j — r sm 6 cos a / - — -r -j — r sm u cos a 



{2nf '""'^ ' J {2nf \p,\ \ p J {27rf \q^\ \ q 



It is easy to see that in the considered symmetrical case when we neglect all contributions generated 
by the tensors g^^, {5i^5jy — 5ij5^i, + 5iy5j^), for the state with the filled in Fermi sphere the angles of 
pairing could be defined by the condition of functional minimum (1181) only for the momenta larger 
than Fermi momentum Pp, see for comparison [5j. The quarks forming the Fermi sphere look like 
the free (non-interacting) ones, see the first term. Now let us calculate the quark chemical potential 
which, by definition, is an energy necessary for adding (removing) one quasi-particle to (from) a 

system jj, = where = 2A'c V I^^ {2n)'^ ~ ^ total number of particles in 



^Just this expression is in the last line of Eq. ([T7]). 
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Pf MeV 

Figure 2: The quark chemical potential as a function of the Fermi momentum for the NJL model. 
The solid line corresponds to the current quark mass m = 5.5 MeV and the dashed one shows the 
behaviour in the chiral limit. 



the volume V. Redefining the chemical potential as /i = f^p . we consider the model with 

correlation function behaving as the S- function in the coordinate space. It is easy to see that 
we come to the Nambu- Jona- Lasinio model (NJL) [B] in this approach. The regularization is 
required to obtain an intelligent result in this model. We adjust the NJL model for the parameter 
set given by Ref. [7], and limit the integration interval over momentum in Eq. (|T71) with the quantity 
IpI < a (A = 631 MeV). Then the functional (ITSl) is written in the following form (inessential terms 
contributing to the constant values are omitted) 



w = Wo + 



A 



Pp {271 



dp 



|n4| (1 — cos6') — G--^ f sin^^ cos^^j / -. — f sin^^' cos 6*' 



(19) 

where Wq = /'^^ ^^^^ contribution coming from the free quarks, m = 5.5 MeV. The 

equation to calculate the equihbrium angle 6 reads as 

(p^ + m^) sin e - Mq{pcos9 + m sin 6*) = , (20) 

where 



Mg = 2G /""t^A (sin^^--cos^^) . (21) 
" (27r)^ |p4| \ P ' ^ ' 



It allows us to obtain a well known selfconsistent gap equation for the dynamical quark mass. 

For the parameters used the dynamical quark mass at zero Fermi momentum is Mq = —335 MeV 
and for the quark condensate 

{N\qq\N) = ^ / dp {p sin 9 - m cos 6) , (22) 

we have (cr|gg|a) = — i (247 MeV)^. The constant characterizing a strength of four-fermion interaction 
was taken as G'A^/(27r^) = 1.34. In Fig. [1] the dynamical quark mass as a function of the Fermi 
momentum is depicted. For comparison the data are presented for current quark mass m = 5.5 MeV 
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(the solid line) and the dashed line corresponds to the calculation in the chiral limit. For the NJL 
model, in particular, the quark chemical potential equals to 



iPfl cosep + Mq 



\P^\ 

K 4 I 



(23) 



When the Fermi momentum reaches zero value the chemical potential quantity coincides with dy- 
namical quark mass /i(0) = \Mg\ — Mg m/\Mq\ = \Mg\ +m. In Fig. [2]the quark chemical potential is 
depicted as a function of the Fermi momentum for the configurations analogous to the ones shown in 
Fig. [H The chemical potential dependence on the Fermi momentum showing up could be interpreted 
as the effect of rapid decrease of the dynamical quark mass with the Fermi momentum increasing. 
Then using Eq. (!20l) the chemical potential may be presented as the following 



Mg Pf 

|Pf |sin^^ 



and taking into account the identity 



P4I' - Mgmf + M2p2 = [p2 ^ _ ^^Y] |P 



41 5 



(24) 



we come to the noteworthy definition of the chemical potential 

^=[Pl + {m-MgfY/K 
Let us remind that for the free fermion gas the chemical potential increases monotonically with the 
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Figure 3: The dynamical quark mass (|Mg|) as a function of the chemical potential. The solid line 
corresponds to the current quark mass m = 5.5 MeV. The dashed one shows the behaviour in the 
chiral limit. 



Fermi momentum growing. The curious feature of the NJL model is the appearance of state almost 
degenerate with vacuum state while the process of filling up the Fermi sphere reaches to the momenta 
close to the dynamical quark mass value (the similar value is characteristic for the momentum of 
quark inside a baryon, see, for example, [3]). This state density with the factor 3 (which expresses the 
relation between baryonic and quark degrees of freedom) absorbed corresponds to a normal nuclear 
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density (n ~ 0.12/fm'^), and chiral condensate could be estimated as | < gg > 



1/3 



100 MeV. In 



the chiral limit the chemical potential near the discussed point is even smaller than the vacuum one. 
The full coincidence of the chemical potentials occurs at the values of current quark mass around 2 
MeV. In fact, Fig. 2 shows that the u quark bond looks stronger than one of the d quark. For clarity 
in Fig. Othe dynamical quark mass, is plotted as the function of chemical potential (Fig. IHfor the 
quark condensate). 



250 




500 

MeV 

Figure 4: The quark condensate (| < gg >^/^ |) as a function of the chemical potential. The solid 
line corresponds to the current quark mass m = 5.5 MeV. The dashed one shows the behaviour in 
the chiral limit. 



The pressure of the quark ensemble 



P = — 



dE 
dV 



dE Pp dE 



'£ + H n 



(25) 



is depicted on the Fig. [5] as a function of the Fermi momentum where n = N/V is the quark density. 
The quark pressure at the values of the Fermi momentum close to the quantity of dynamical quark 
mass is approximately degenerate with the vacuum pressure (slightly lower than the vacuum one). 
The vacuum density is of order 40 — 50 MeV/fm^ and corresponds well to the value extracted from 
the bag models, see, for example, Ref. [3|. Actually, all the NJL results could be received in the 
mean field approximation because the trigonometrical terms in the mean energy definition Eq.( fT9l) 
can be rewritten (using Eq. (!24l) again) as the functions of dynamical and current quark masses in 
the following form 



p sin 6 — m cos 6 

\Pi\ 

\p/i\ cos 9 



Mg-m 



[p^ + (Mq - m 
+ m{m — Mq) 
[p2 + (M,-m)2]V' ■ 



,211/2 



In order to trace back the dependence of all results on the formfactor form we consider the model 
(in a sense, opposite to the NJL model) with the formfactor behaving as 5-function in momentum 
space I{p) = (27r)'^ ^(p)- This limit is an analogue of the Keldysh model which is well known in the 
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Figure 5: The pressure of the quark ensemble as function of the Fermi momentum. The sohd hne 
for the current quark mass m = 5.5 MeV. The dashed one — in chiral hmit. 



physics of condensed matter [8] and the mean energy functional f|T8l) develops the following form 



w 



Pf dp 



ml + 



dp 



\Pa\ (1 



cos( 



G 



dp 



i2nf Jp, {2nf ^^^^^ ^ Jp, (27r)3 \p,\^ 

The chemical potential is defined in this approach as 

[Pp sin Op — m cos OpY 



sin 6* — — cos 

K P / 



(26) 



/i 



IP4 I cosOp 



G- 



I P-F|2 
K 4 I 



(27) 



It results from this definition that at low Fermi momenta the chemical potential goes to fi —>■ m + G. 
In Ref. the constant G was taken of order of dynamical quark mass in NJL. With Fermi momentum 
increasing the chemical potential remains approximately constant and starts to increase when the 
Fermi momentum exceeds the G value. The vacuum density (Pp = 0) turns out to be singular and 
for the pressure difference we obtain 



P-P(O) 



2 N, 
27r2 



dp p cos6'|p4| —G dp 



Pf 



p 



b4p 



sin 6 cos ( 

V P 



:28) 



This function is slowly monotonically growing till the Fermi momentum value of G. 

In summary we would like to emphasize our main result which, as we believe, is quite transparent. 
Our estimate of the effects responding to the process of filling up the Fermi sphere demonstrates 
the parallels between this quasi-particle picture and the conventional bag model but with one new 
essential element. It is just the presence of instability region dP/dPp < 0. Then the states ([6]) could 
be considered as natural 'building' material for baryons. In principle, one of the ways to construct 
the corresponding bound state could follow the Walecka model ideas [U| with utilizing the information 
on scalar and vector field behaviour and respective constants of their interaction with quarks which 
have been gained up to now [5]. 

This work was supported by the INTAS Grant 04-84-398 and the Grant of National Academy of 
Sciences of Ukraine. 
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